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ABSTRACT 

This article presents two new automation tools that 
obtain stresses and strains (sheai and peel) in adhe- 
sively bonded joints. For a given adhesively bonded 
joint Finite Element model, in which the adhesive is 
characterised using springs, these automation tools 
read the corresponding input and output files, use the 
spring forces and deformations to obtain the adhesive 
stresses and strains, sort the stresses and strains in 
descending order, and generate plot files for 3D visu- 
alisation of the stress and strain fields. Grids (nodes) 
and elements can be numbered in any order that is 
convenient for the user. Using the automation tools, 
trade-off studies, which are needed for design of adhe- 
sively bonded joints, can be performed very quickly. 

Keywords: Finite Element Analysis, Adhesively 

Bonded Joints, Automation, Stiffness, Shear/Peel 
Stress/Strain 

1 NOMENCLATURE 

All of the symbols used in this article are defined in 
Table 1. 


Table 1: Definition of symbols. 


Symbol 

Definition 

A 

area of plate elements 
used to model adherends 

E 

adhesive elasticity modulus 

fx,l 

X spring force at grid i 

Am 

Y spring force at grid i 

fz,i 

Z spring force at grid i 

G 

adhesive shear modulus 

i 

grid number 


Table 1: Definition of symbols (continued). 


Symbol 

Definition 

K p 

peel spring stiffness 

K s 

shear spring stiffness 

P 

peel force 

V 

shear force 

A 

peel deflection 

8 

shear deformation 


X spring deformation at grid i 


Y spring deformation at grid i 

A,* 

Z spring deformation at grid i 

e 

adhesive peel strain 

€i 

peel strain at grid i 

7 

adhesive shear strain 

7i 

shear strain at grid i 

7] 

thickness of adhesive 

V 

adhesive Poisson's ratio 

a 

adhesive peel stress 

Oi 

peel stress at grid i 

T 

adhesive shear stress 

Ti 

shear stress at grid i 


2 INTRODUCTION 

L.J. Hart-Smith has performed a lot of pioneering 
work in analysis and design of adhesively bonded 
joints (e.g., [1], [2]). 

Many others have used Finite Element methods 
and springs (linear or non-linear) to model adhesives 
in bonded joints (e.g., [3]). The Finite Element model 
of a typical adhesively bonded joint is shown in Fig- 
ure 1. If forces and/or moments are applied to such 
a joint, spring forces and deformations can be used 
to determine the stresses and strains in the adhesive. 




This document describes two FORTRAN programs 



Figure 1: Finite Element model of a bonded joint. 


which automate the process of obtaining stresses and 
strains in adhesively bonded joints. Given an adhe- 
sively bonded joint model, these programs read the 
corresponding NASTRAN (NAS A Structural Analy- 
sis, a Finite Element software package widely used in 
the Aerospace industry) input and output files, use 
the spring forces or deformations to obtain the adhe- 
sive stresses or strains at the mid-plane grid points, 
sort the stresses and strains in descending order, and 
generate Mathematica (a popular symbolic mathe- 
matics and graphics software package) plot files for 
3D visualisation of the stress and strain fields. 

Note that the grids and springs can be num- 
bered in any order. The programs determine the 
pairs of grids which are initially -coincident, and the 
spring triplets which correspond to such grid pairs. 
In the earlier attempts at NASA GSFC to model ad- 
hesives in bonded joints, the analyst had to number 
the grids and/or springs in specific orders so that the 
XYZ spring triplets acting between pairs of initially- 
coincident adhesive grids could be identified. Such 
ordering of grids and/or springs can be very time- 
consuming. 

The spring stiffness values are determined from the 
adhesive stress-strain curve and the bonded joint ge- 
ometry (see Section 4). 


3 THE SPRING MODELING 
PROCESS 

The modeling process involves the following steps: 

1. A fine mesh of grids is created in the mid-plane 
of the bonded joint. 

2. For every grid in the mid-plane mesh, an 
i nitially-coincident grid is created. 

3. Three springs are placed between each pair of 
initially-coincident grids. These springs act in 
the X, Y, and Z directions; where the X and 
Y axes are parallel to the overlap plane of the 
bonded joint, and the Z axis is defined by the 
right-hand- rule. 

4. Rigid elements are used to connect the mid- 
plane grids to the corner grids of the plate el- 
ements which represent the adherends (see Fig- 
ure 2). 
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Figure 2: Depiction of shear and peel springs. 

4 DETERMINATION OF THE 
SPRING STIFFNESS VALUES 

4.1 Stiffness of Shear Springs 

To obtain the shear spring stiffness values, we can 
use the following general relationships: 

y 

t = Gy, t = — , 5 = 7 V l 1 ) 




Symbols in equation (1) are defined in Table 1. 
Therefore, 



We also know that: 

V = K s S (3) 


practice), then the relationships specified in equa- 
tion (5) apply, and must be used in calculating the 
peel spring stiffness values. 

If the adhesive elasticity module is not published 
in the manufacturer’s literature, it can be calculated 
from its Poisson’s ratio and its shear modulus using 
the following equation: 


where K s is defined in Table 1. Combining equa- 
tions (2) and (3), we get 


K s 


GA 


(4) 


E = 2G{l + v) - (10) 

Loss and Keyward [3] discuss variations of equa- 
tion (10) for various stress states. 


Equation (4) gives the shear spiing stiffness values 
as a function of the adhesive shear modulus and the 
model geometry, and can be used in assigning stiff- 
ness values to the shear springs of an adhesive Finite 
Element Model. Note that if quadrilateral shell ele- 
ments are used to model the adherends (a common 
practice), then: 

A = Asternal = 2 A e d ge = 4 A corneT (5) 

where Ai n t erna h A e dgei & n d A corni[T represent the ele- 
ment areas to be used in equation (4) for calculating 
the shear spring stiffness values. 


4.2 Stiffness of Peel Springs 


To obtain the peel spring stiffness values, we can use 
the following general relationships: 


P 

a — Ee, a = — , 


A 

rj 


(6) 


5 DESCRIPTION OF THE AD- 
HESIVE STRESS PROGRAM 

Given a bonded joint model, the program executes 
the following operations. 

1. Prompts the user for the name of NASTRAN 
input and output files. 

2. Locates and stores all of the spring ID’s, the 
corresponding grid ID’s, and the corresponding 
grid coordinates in the NASTRAN input file. 

3. Identifies the pairs of grids which are coincident 
in the unloaded model. 

4. Identifies the spring triplets (X, Y, and Z) for 
the pairs of initially-coincident grids. 

5. Locates and stores the spring forces in the cor- 
responding NASTRAN output file. 


Symbols in equation (6) are defined in Table 1. 
Therefore, 

- = E- (7) 

A t 7 


We also know that: 


P = K p A 


( 8 ) 


where K p is defined in Table 1 Combining equa- 
tions (7) and (8), we get 


K p = 


EA 


(9) 


6. Assigns the stored spring forces to the appropri- 
ate springs in the triplets identified in step 4. 

7. Determines adhesive shear and peel stresses at 
the mid-plane grids using the following equa- 
tions. 


n = 


//■2 , f 2 . 

V J x,i 

A 


<7, = 


fz,i 


(11) 

( 12 ) 


The symbols in equations (11) and (12) are de- 
fined in Table 1. 


Equation (9) gives the peel spring stiffness values as 
a function of the adhesive elasticity modulus and the 
model geometry, and can be used in assigning stiff- 
ness values to the peel springs of an adhesive Finite 
Element Model. Note that if quadrilateral shell ele- 
ments are used to model the adherends (a common 


8. Sorts the shear and peel stresses in descending 
order and writes them to an output file. 

9. Writes the coordinates of the mid-plane grid 
points and their corresponding shear and peel 
stresses to plot files. 



The plot files mentioned in step 9 are used to gen- 
erate 3D representations of the si ear and peel stress 
fields using Mathematica. Figure 3 shows the 3D 
shear stress field for the joint represented in Fig- 
ure 1. To efficiently execute the operations spec- 


Tension load of 16< 00 lb. 
is applied in the / direction. 



2. Calculates and stores the adhesive spring defor- 
mations from the grid point displacements ob- 
tained in the previous step. 

3. Assigns the stored spring deformations to the 
appropriate elements in the spring triplets corre- 
sponding to the pairs of initially-coincident grids 
(see step 4 of the previous section). 

4. Calculates adhesive shear and normal strains at 
the mid-plane grids using the following equa- 



The symbols in equations (13) and (14) are de- 
fined in Table 1 . 

5. Sorts the shear and normal strains in descending 
order and writes them to an output file. 


Figure 3: Shear stress in a bonded joint. 

ified in steps 3, 4, and 6; the program invokes var- 
ious sorting subroutines. These sorting operations 
and other details of the program can be found in 
the source code (adhsv_strs.f) which is located at the 
http://femci.gsfc.nasa.gov/adhes ve site. 


6. Writes the coordinates of the mid-plane grid 
points and their corresponding shear and nor- 
mal strains to plot files. 

The shear and normal strain plot files are read 
into Mathematica to generate 3D representation of 
the strain fields. Figure 4 shows the 3D shear 
strain field for the bonded joint represented in 
Figure 1. The program invokes various sort- 


6 DESCRIPTION OF THE AD- 
HESIVE STRAIN PROGRAM 

Equations (11) and (12) show that the adhesive shear 
and peel stresses are dependent on the area of the 
plate elements used to represent the adherends. For 
complicated joints, it may not be possible to have a 
uniform area for such plate elements. Therefore, as 
an alternative, a program was written to determine 
shear and peel strains which are independent of the 
plate element area, and dependent on the adhesive 
thickness which is normally uniform. 

For a given bonded joint model, the program starts 
by executing steps 1 through 4 that the Adhesive 
Stress program goes through (see the previous sec- 
tion). Next, the program performs the following op- 
erations. 

1. Locates and stores the mid- plane grid point dis- 
placements in the N ASTRA N output file. 



Figure 4: Shear strain in a bonded joint. 

ing subroutines to efficiently execute some of the 
steps specified above. These sorting operations and 
other details of the program can be found in the 
source code (adhsv_strn.f) which is located at the 
http://femci.gsfc.nasa.gov/adhesive site. 


7 SAMPLE OUTPUT FILES 

7.1 Adhesive Stress Program Output 

A short portion of a sample output file for the adhe- 
sive stress program is shown below. 

No. of springs = 1683. 

No. of points = 561 

Springs: 4002 5683 5684 => Tau= 17050.4878, 

Sigma= 5020.4877 
Springs: 4005 5685 5686 => Tau= 21762.8052, 

£ igma= -3518.7251 

Springs: 4008 5687 5688 = > lau= 14779.8012, 

Sigma= -6167.9155 


********* SORTED SHEAR STRESSES ********* 
Springs: 4194 5811 5812 = > lau = 38807.8485 
Springs: 5580 6735 6736 => Tau= 38807.8485 
Springs: 5490 6675 6676 => 1au = 38713.2117 


********* SORTED PEEL STRESSES ********* 
Springs: 4200 5815 5816 => Sigma= 15544.5325 
Springs: 5388 6607 6608 => Sigma= 15544.5325 
Springs: 5289 6541 6542 => Sigma= 15004.5107 


7.2 Adhesive Strain Program Output 

A short portion of a sample output file for the adhe- 
sive strain program is shown below. 

No. of spring grids = 561 
Grids: 8001 9001 => Gairana= 0.3811, 

Epsilon= -0.1195 

Grids: 8002 9002 => Gamma= 0.2440, 

Epsilon= 0.0419 

Grids: 8003 9003 => Gamma= 0.1677, 

Epsilon= 0.0734 


***** SORTED SHEAR STRAINS ***** 
Grids: 8033 9033 => Gamma= 0.3848 
Grids: 8561 9561 => Ganuna= 0.3848 
Grids: 8001 9001 => Gamma= 0.3811 


******* SORTED PEEL STRAINS ******* 
Grids: 8099 9099 => Epsilon= 0.1797 
Grids: 8495 9495 => Epsilon= 0.1797 
Grids: 8462 9462 => Epsilon= 0.1786 


8 SUMMARY 

In this article, two new automation tools are intro- 
duced and described. 

The automation tools are used to quickly deter- 
mine shear and peel stresses and strains in adhesively 
bonded joints. As a result, trade-off studies for de- 
sign of such joints can be performed rapidly. 

Given an adhesively bonded joint NASTRAN 
model, in which the adhesive is characterised using 
springs, the two automation tools read the corre- 
sponding input and output files, use the spring forces 
and deformations to obtain the adhesive stresses and 
strains, sort the stresses and strains in descending 
order, and generate plot files for 3D visualisation of 
the stress and strain fields. 

The automation tools can handle any scheme used 
by the user to number the model elements and grids 
(nodes). 

The results of both programs have been verified 
by hand-calculating stresses and strains at various 
adhesive grid points. 
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